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Методами моделирования представительных объёмов, эффективных модулей и конечных элементов опреде- 
ляется полный набор материальных констант пористых пьезокерамических материалов с учётом гипотезы о 
неоднородной поляризации в окрестности пор. Для численного определения эффективных модулей в конеч- 
но-элементном пакете АМ5У$ решаются наборы статических пьезоэлектрических задач для представитель- 
ных объёмов с граничными условиями, обеспечивающими постоянные значения полей деформаций, напря- 
жений, электрического поля и электрической индукции однородного тела. Расчёты проводятся с помощью 
специально разработанных компьютерных программ, написанных на макроязыке АРОЁЕ АМ5У$. Исследуется 
влияние количества ячеек и различных моделей представительных объёмов на значения важнейших эффек- 
тивных модулей для ряда типов пористой пьезокерамики. Обсуждается корреляция результатов вычисли- 
тельных и натурных экспериментов. Показано, что характер зависимостей эффективных модулей от пори- 
стости слабо зависит от типа пористых пьезоэлектрических материалов. 
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Введение. В настоящее время пористая пьезокерамика, являющаяся сравнительно новым клас- 
сом пьезоэлектрических композиционных материалов, находит широкое применение в ультразву- 
ковых преобразователях, гидрофонах, датчиках давления и других пьезоэлектрических устрой- 
ствах. Основная выгода при использовании пьезокомпозитов заключается в возможности суще- 
ственного улучшения основных рабочих параметров пьезоустройств, особенно при их нагрузке на 
акустические среды. Так, экспериментальные исследования показывают [1, 2], что пористая пье- 
зокерамика обладает высокой объёмной пьезочувствительностью в широкой полосе частот и бо- 
лее низким импедансом по сравнению с плотной пьезокерамикой. При этом указанные эффекты 
усиливаются при повышении степени пористости керамики. В то же время важнейшие толщинные 
характеристики (пьезомодуль @.,, коэффициенты электромеханической связи К,, К) для ряда 
пористых пьезокерамик практически не зависят от пористости, а соответствующие продольные 
величины (4.,, К,, К,, ) быстро убывают с ростом пористости [1—3]. Между тем, имеющиеся экс- 


периментальные данные по свойствам пористой пьезокерамики достаточно противоречивы, и за- 
висимости пьезомодулей, например в [4—6], отличаются от полученных в [1—3]. Всё это опреде- 
ляет актуальность решения задач по теоретическому определению и предсказанию свойств пье- 
зокомпозиционных материалов, и, в частности, пористых пьезокерамик. Обзоры по разнообраз- 
ным исследованиям пористых пьезокомпозитов можно найти в [1, 7] и в многочисленных статьях 
А. Н. Рыбянца. В представленной работе рассматривается один из методов механики компози- 
тов — метод эффективных модулей [3, 7—9]. Уточнённое моделирование базируется на исполь- 
зовании специальных моделей представительных объёмов, метода конечных элементов (МКЭ) и 
учёта неоднородности поляризации пьезокерамики вблизи пор. Результаты вычислительных экс- 





1 Работа выполнена при частичной поддержке АВЦП «Развитие научного потенциала высшей школы» Минобрнауки (про- 
ект 1.6052.2011) и РФФИ. 
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периментов, рассчитанных по представленной технологии моделирования, призваны ответить на 
вопрос о том, существенно ли меняются теоретические зависимости пьезомодулей и других эф- 
фективных характеристик от пористости для пьезокерамик различной сегнетожёсткости. 

Базовый вариант метода эффективных модулей. Пусть @ — представительный объём не- 
однородного пьезокерамического материала кристаллографического класса 6бтт, Г=0® — его 
граница. Рассмотрим следующую систему дифференциальных уравнений статической теории 
электроупругости в ©: 


Е(У).Т=0, У.0=0 
Т=с.5-е’.Е, О-е.5+=.Е 
$=ЦУ)-ц, Е=-Уф 
(У) (1) 
обсбаа д, 
Е(У=0 ба бд, У=:6, 
обаоаао д, 
с граничными условиями 
и=Е(х).$,, Ф=-х. Е, хеГ=2О. (2) 


В (1)—(2) для удобства дальнейшего применения метода конечных элементов использованы век- 
торно-матричные обозначения, стандартные для технических применений теории электроупруго- 


сти (пьезоэлектричества): Т ={0,,,0,„,0..,0,.,0,.,01›! — псевдо-вектор напряжений о,; В — 


вектор электрической индукции; $ = Ее, 2Еь | — псевдо-вектор деформаций =, ; 
Е — вектор напряжённости электрического поля; и=и(х) — вектор-функция перемещений; 
Ф=Ф(х) — функция электрического потенциала; с“ = с" (х) — матрица упругих жёсткостей, 
измеренных при постоянном электрическом поле, размерности 6х6; е=е(х) — матрица пьезо- 
модулей, размерности 3хб; =° = =° (х) — матрица диэлектрических проницаемостей, измерен- 
ных при постоянных деформациях, размерности 3х3; х — вектор пространственных координат; 
(...)' — операция транспонирования для матриц и векторов. 

Особенностью краевой задачи (1)—(2) являются граничные условия (2), в которых пред- 
полагается, ЧТО! $; = {оао Вена р Во = Е.Е» Е ГДЕ ВОВЕ ур == ПОСТО- 
янные (не зависящие от х ) величины. Для такой задачи, как известно [6], ($) =$5., (Е) =Е, , где 


в угловых скобках обозначаются осреднённые по объёму характеристики: Са) — 5 |. (-..)а® . Кро- 


ме того, будут выполняться и важные энергетические соотношения: (5 т) = (5") . (т) = $, . (т) ; 


(Е'.0) = (Е°). (2) = Е, - (0). 

Указанные свойства позволяют обосновать метод нахождения эффективных модулей 
С", е*", =5*", известный как метод эффективных модулей в механике композитов [3, 7—9]. 
Поставим в соответствие неоднородному пьезоэлектрическому композиту в области @ среду 


сравнения — макрооднородный пьезоэлектрический материал с некоторыми неизвестными пока 
модулями с““", е*", =5*". Будем рассматривать для среды сравнения в © ту же задачу (1)—(2), 


но с постоянными модулями с““", е*", =5“". Легко проверить, что решением такой задачи в об- 
ласти © будут являться следующие механические и электрические поля: 
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и=и, =Ё (х).5,, Фф=ф=-х-Е,, $=5$,, Е=Е,, (3) 
Т=Т=с““".$, -е*".Е, б=б, =е* .$, + =5".Е, (4) 
причём векторы механических напряжений Т =Т,, деформаций $ =$,, электрической индукции 


О =0, и напряжённости электрического поля Е =Е, не зависят от пространственных координат 
х. Если предположить, что для неоднородного пьезокомпозита и для среды сравнения равны 
средние напряжения и средняя электрическая индукция: 

(т)= (то), (2) = (В), (5) 
то это обеспечит и равенства средних деформаций, векторов напряжённости электрического поля 
и потенциальных электромеханических энергий П = (5 -Т+Е“ :0) | 


Соотношения (4), (5) позволяют получить расчётные формулы для определения эффек- 
тивных модулей. Именно, как несложно заметить из (4), (5), для этого достаточно в задаче (1)— 
(2) для неоднородного композита положить в краевых условиях (2) лишь одну из компонент $,., 
Е, (а=1.....6;/=1,2,3) отличной от нуля. Тогда, предполагая, что эффективная среда сравне- 
ния также является пьезокерамическим материалом класса бт, можно выделить пять краевых 
задач (1)—(2), различающихся только ненулевыми компонентами в граничных условиях (2), кото- 
рые позволяют в совокупности найти полный набор эффективных констант пьезокерамического 
композита: 

1. $, =&, 5, =0, В=1, Е; =0 => 





с" = (0,6, 1=1,2,3, ей = (В; в. ; (6) 
П. $; =&, 5, =0, В*З, Е, =0 > 
ты 7 (0.3 ) /Е ' а г. (2; Е ' о = (о: )/Ео я (о„› ) /Е р (7) 
Ш. $, =25,, $3 =0 В, Е == 
бт = (0 Ека =, ) (25, (8) 
М. $, =0, Е, = Ее, 
в =-(о..)/Е. ' ИА = (2 ИЕ р (9) 





М $, =0, Е, = Ее, 


ез =- (о, )/Е, =-(о„,)/Е, ' ез =- (0. )/Е, й Е = (р, Ес у (10) 


Итак, чтобы определить десять независимых эффективных констант пористой пьезокерамики 


Е ей Е её Е ей Е ей Е ей у у у 5 ей 5 ей 
(Си, С”, Св’, Св, Си ‚ел, е, е5, =”, & ), достаточно решить пять статических 


задач электроупругости (1)—(2) с различными $, и Е, в краевых условиях (2) и использовать 
формулы (6)—(10). 
Моделирование представительных объёмов. Приведённые формулы расчёта эффективных 
модулей предполагают решение соответствующих краевых задач электроупругости в областях ©, 
которые должны являться представительными объёмами пористых пьезокомпозитов. В качестве 
представительных объёмов в идеале следует брать области, достаточно большие по сравнению с 
размерами неоднородностей (т. е. пор), но малые по сравнению с расстояниями, на которых су- 
щественно меняются медленные переменные. 

Как известно [1], пористый пьезокерамический материал является бинарным композитом, 
состоящим из пьезоактивной и пассивной (непьезоактивной) фаз. Пьезоактивная фаза образует 
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трёхмерный пьезокерамический каркас, в котором пассивная фаза существует в виде изолиро- 
ванных или связанных друг с другом воздушных пор. Первый случай по классификации 
Р. Э. Ньюнема отвечает связности 3—0, а второй — связности 3—3 (закрытая и открытая пори- 
стость соответственно). 

При небольшом проценте пористости достаточно простой, но адекватной наблюдаемой на 
практике микроструктуре пористого материала, представляется модель пьезокерамической куби- 
ческой решётки, состоящей из одинаковых ячеек-кубиков, часть из которых случайным образом 
объявляется порами. Отметим, что подобная модель рассматривалась в [3]. Однако такая модель 
при большом числе пор может потерять связность каркаса, а также не поддерживает структуру 
связности композита (3—0 или 3—3). 

Для построения связанных структур в кубической решётке можно использовать алгоритмы 
теории перколяции, позволяющие получать протекающие кластеры. При малой пористости можно 
строить кластеры из пор, а при большой пористости — кластеры из пьезокерамического материа- 
ла. Один из таких методов — метод ограниченной диффузией агрегации Виттена — Сандера [10] 
был программно реализован В. В. Ремизовым и наряду с другими аналогичными методами анали- 
зировался в [11] применительно к рассматриваемым задачам. 

Конечно-элементное решение задач для представительных объёмов. Для решения задач 
электроупругости (1)—(2) для неоднородного пьезоэлектрического материала в представитель- 
ном объёме © можно перейти к их слабым постановкам и использовать классическую технику ко- 
нечно-элементных аппроксимаций. Пусть ©, — область, занимаемая соответствующей конечно- 
элементной сеткой, 9, < ©, 9, =Ц, 9“, где 9%“ — отдельный конечный элемент с номером К. 
Для статических задач можно искать приближённое решение {и, =им, Ф, =Ф} на конечно- 
элементной сетке ©, =, 9“ в форме 

и, (х,Ё) = М, (х)- Ц, ф, (х,Ё) =№ (х).Ф, (11) 


где № — матрица базисных функций для перемещений; М, — вектор-строка базисных функций 


для электрического потенциала; М, Ф — векторы узловых перемещений и электрического 
потенциала. 

В соответствии с обычной конечно-элементной техникой, аппроксимируем континуальную 
слабую постановку задачи электроупругости в конечномерных пространствах, связанных с базис- 
ными функциями №, М, . Подставляя (11) и аналогичные представления для проекционных 
функций в слабую постановку задачи электроупругости для ©®,, получим следующую конечно- 
элементную систему: 

в (12) 


-К.,-9+К„-Ф-Е, 
Здесь К„=»У“К%, К, =У“”КХ, К» =>"К* — глобальные конечно-элементные матрицы, 


а 
получаемые из элементных матриц в результате процесса ансамблирования Г» РР. ве 


торы, определяемые внешними воздействиями. 
Элементные матрицы имеют вид: 


кх = [|.5=.с'.5°09, К = [.52.е”.5°49, (13) 


К = [1 $5 -=°. 5500, $8 = МУ) №”, $5 =У№", (14) 
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где №”, №" — матрица и вектор-строка аппроксимирующих базисных функций, соответственно 


определённых на отдельных конечных элементах. 
Отметим, что МКЭ является сейчас стандартным методом решения задач пьезоэлектриче- 
ского приборостроения и с успехом используется как в практике расчётов реальных пьезоэлек- 
трических устройств, так и на микро- и наноуровнях исследований в электроупругости [12]. 
Моделирование неоднородной поляризации. При рассмотрении неоднородного пористого 
пьезоэлектрического материала логично предположить, что поляризация пьезокерамики может 
быть неоднородной и наследовать структуру поля поляризации. В связи с этим, для учёта частич- 
ной поляризации пьезокерамики в окрестности пор, на предварительном этапе можно промоде- 
лировать процесс поляризации пористой керамики. Для этого решим задачу квазиэлектростатики 
с помощью МКЭ для пористого диэлектрика в области © , сгенерированной для представительно- 
го объёма по одному из двух методов, описанных в предыдущем разделе. При этом каждую из 
кубических ячеек можно считать кубическим конечным элементом (КЭ) со своими диэлектриче- 
скими свойствами, различными для КЭ, заполненных материалом, и для пористых КЭ. 
Тогда для неоднородного куба © =[0,/ |х| 0,1 |х| 0, | в декартовой системе координат 
Ох, хх, , будем иметь следующую краевую задачу: 
У.0=0, О==.Е, Е=-Уф, ХЕО, (15) 
ФЕИ,, хеЕГы, Л=12; п.9=0, хеГ,, (16) 
где [ — линейный размер куба; Г=Ц,Г.; УГ. ; Го; 


пропорциональная значению поля поляриза- 


— электроды х, =0 и х, =[, для которых 
задаётся разность потенциалов И, —И, = к Е „, 
ции Е, с коэффициентом пропорциональности к; == =(х) — матрица диэлектрических про- 


ницаемостей неполяризованной керамики с порами. 
Заметим, что система МКЭ для задачи (15), (16) получается фактически как частный слу- 





чай из (12)—(14), если положить с" =0, е=0, =° ==, ине использовать степени свободы Ц. 
После решения задачи (15), (16) находились значения векторов поляризации 





Р* =0“ — = Е* в центральной точке каждого КЭ с номером А, не являющегося порой. Здесь и 
далее через =, обозначена диэлектрическая проницаемость вакуума, =, = 8,85.10-? Ф/м. С дан- 
ными КЭ ассоциировались свои элементные системы координат Ох“х*х=, для которых оси 
Ох* выбирались такими, чтобы их направления совпадали с направлениями векторов поляриза- 
ции Р“. 

Две другие оси определялись по следующему алгоритму (рис. 1). Пусть Р* — вектор про- 
екции вектора поляризации Р“ на плоскость Ох,х, исходной декартовой системы координат, 


т.е. Р® =Р“ - Р*е,, где е, — орт оси Ох,. Первым поворотом относительно оси Ох, в плоско- 
сти Ох.х, на угол 0% повернём оси Ох,, Ох, так, чтобы новая ось ОХ” совпала по направле- 
нию с вектором Р*. Получив промежуточную систему координат Ох*х”х, , повернём её вторым 
поворотом относительно оси ОХ“ в плоскости Ох,Х“ на угол 05 так, чтобы ось Ох“ совпала 
р*^ 


по направлению с вектором поляризации Р“, то есть со50% = Р* / ‚ Охх =0Х“ . В резуль- 








тате получим требуемую элементную систему координат Ох*х*х=. 
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Рис. 1. Элементная система координат отдельного конечного элемента 


Для работоспособности указанного алгоритма при малых величинах модуля вектора Р®, 


то есть в случае, когда направление вектора поляризации Р“ близко к направлению оси Ох,, 


алгоритм модифицировался следующим образом: со$0% = Р*° | Р> |, если |Р®|> 4%; с0$0*% =1, 














если Р®' |< 4, где 4% =0,001|Р* 

На втором этапе моделирования КЭ электростатики модифицируются в элементы с воз- 
можностями пьезоэлектрического анализа. Новым элементам присваиваются материальные свой- 
ства двух типов: поляризованной пьезокерамики для КЭ матрицы и пренебрежимо малые упругие 
и пьезоэлектрические модули и диэлектрические проницаемости, равные диэлектрической про- 


ницаемости вакуума, для пор. С поляризованными КЭ пьезокерамической матрицы связываются 














элементные системы координат Ох“х*х= , определённые по векторам поляризации Р“. Далее 


для определения эффективных модулей решается задача электроупругости (1)—(2) для неодно- 
родного представительного объёма по методам, описанным выше. Подчеркнём, что при учёте не- 
однородной поляризации задача (1)—(2) решается для неоднородной структуры, причём каждый 
КЭ поляризованной пьезокерамики имеет свои модули с“, е*, =5* , получаемые по известным 
формулам пересчёта тензорных коэффициентов при переходе от кристаллографической декарто- 
вой системы координат Ох,х›х, в элементные системы координат Ох“х*х*. 

Если не учитывать неоднородность поляризации, то задача (15), (16) не используется, и в 
задаче (1)—(2) все КЭ имеют или свойства пьезокерамического материала класса бт, поляри- 
зованного вдоль оси Ох, , или свойства пор. 

Результаты расчётов и выводы. По описанной выше методике были проведены численные 
эксперименты по определению эффективных модулей для пористых пьезокерамических материа- 
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лов различной сегнетожёсткости: сегнетожёсткой пьезокерамики Р7Т-8, пьезокерамики РЕТ-4 
средней сегнетожёсткости и сегнетомягкой пьезокерамики РЕТ-5Н. Входными данными для расчё- 
тов принимались значения модулей пьезокерамических материалов при нулевой пористости, при- 
ведённые в табл. 1. 
Таблица 1 
Модули пьезокерамических материалов различной сегнетожёсткости 
при нулевой пористости 



































Р7Т-8 Р7Т-4 РЕТ-5Н 
СЕ ‚Га 146,8 138,9 1272 
СЕ ‚Га 81,08 77,83 80,2 

© та 81,05 74,28 84,67 
сЕ ‚ГПа 131,7 115,4 117.4 
се ‚ГПа 31,35 25,64 22,98 
е; , Кл/м” 3,875 -5,202 —6,623 
езз , Кл/м? 13,91 15,08 23,24 
е,; ‚ Кл/м? 10,34 12,72 17,03 
5 904 762 1704,4 
25. /&. 561 663 1433,6 




















Конечно-элементные модели представительных объёмов пористых пьезокерамик состояли 
из кубических решёток геометрически одинаковых ячеек-кубиков, которые являлись гексаэдраль- 
ными восьмиузловыми КЭ с пьезоэлектрическими свойствами. В зависимости от заданной величи- 
ны объёмной пористости р часть КЭ объявлялась порами, причём рассматривались два метода 
генерирования пор: простой случайный метод и метод Виттена — Сандера. 

Кубические решётки брались размера лхлхл , где п — число ячеек вдоль одной из коор- 
динатных осей. Вычислительные эксперименты с числами п= 10, 20, 30, 40 показали, что значе- 
ния эффективных модулей стабилизируются при п=20. Это число п=20 и было взято за базо- 


вый размер кубической решётки при дальнейших расчётах. 


После вычисления полных наборов эффективных модулей с„“", е\', &’“" по известным 


формулам находились и другие характеристики пьезокерамических материалов: коэффициенты 
электромеханических связей, упругие податливости, пьезомодули иных форм определяющих со- 


отношений и т. п. Проанализируем поведение пьезомодулей &4“" и а=" в зависимости от пористо- 
сти р. Чтобы сравнить, как изменяются эти пьезомодули для пьезокерамических композитов раз- 
ЛИЧНОЙ сегнетожёсткости, будем рассматривать их относительные изменения 
г(9,,) =а57 (р)/а., ‚где 7=1 или 7=3; 45 (р) — значения эффективных пьезомодулей при 
пористости р; а;, — значения пьезомодулей плотной пьезокерамики, то есть при р =0. 

На рис. 2 приведены зависимости от пористости р относительных величин пьезомодулей 
г (а, ) и г (аз } для сегнетомягкой пьезокерамики РЕТ-5Н в предположении однородной поляри- 


зации пьезокерамической матрицы для двух методов генерации представительных объёмов: про- 
стого случайного метода и метода Виттена — Сандера. Как проиллюстрировано на рис. 2, а, при 
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методе Виттена — Сандера с ростом пористости наблюдается большее убывание пьезомодуля 
ал' (р) ‚ чем при простом случайном методе. 


4 
0.95 


0.9 
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Рис. 2. Зависимости относительных пьезомодулей (1) (а) и Каз) (6) от пористости для сегнетомягкой пьезокерамики 
РЕТ-5Н при однородной поляризации матрицы для различных методов генерации представительных объёмов 


В то же время величины пьезомодуля 45’ (р) практически не зависят от пористости для 


обоих методов (рис. 2, 6). Аналогичные тенденции наблюдались и для пьезокерамики Р7Т-4 сред- 
ней сегнетожёсткости и для сегнетожёсткой пьезокерамики Р7Т-8. 

Таким образом, способы генерации представительных объёмов могут оказывать достаточ- 
но существенное влияние на значения эффективных модулей пористых пьезокомпозитов различ- 


ной сегнетожёсткости. Поскольку большее убывание пьезомодулей аз’ (р) пористых пьезокера- 
мик было получено и в результате обработки ряда натурных экспериментов [1—6], то следует 


признать метод Виттена-Сандера генерации представительных объёмов более подходящим в 
сравнении с простым случайным методом. 


На рис. 3 и рис. 4 отражены зависимости относительных величин пьезомодулей г (а.,) и 
г(а,.) от пористости р для пьезокерамических материалов различной сегнетожёсткости в пред- 


положениях однородной (а) и неоднородной (6) поляризации пьезокерамической матрицы при 
генерации представительных объёмов по методу Виттена — Сандера. 


ТЕ, 





$ $ Сетнетожесткая керамика 
&4& Керамика средней сегнетожесткости 
ФеФ Сегнетомягкая керамика 


$ $ Сетнетожесткая керамика , 
444 Керамика средней сегнетожесткости | *. 
ФФ Сегнетомягкая керамика 
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Рис. 3. Зависимости относительных пьезомодулей (1) от пористости для пьезокерамик различной сегнетожёсткости 
при однородной (а) и неоднородной (6) поляризациях матрицы при генерации представительных объёмов по методу 
Виттена — Сандера 
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Рис. 4. Зависимости относительных пьезомодулей /(%&з) от пористости для пьезокерамик различной сегнетожёсткости 
при однородной (а) и неоднородной (6) поляризациях матрицы при генерации представительных объёмов по методу 
Виттена — Сандера 


Как видно на рис. 3 и 4, для более сегнетомягких керамик (т. е. для пьезокерамик с боль- 
шей подвижностью доменных стенок) пьезомодули 4з' (р) и а; (р) несколько сильнее зависят 


от пористости. Однако основные тенденции поведения пьезомодулей от пористости сохраняются, 
а сами относительные величины пьезомодулей достаточно близки для всех рассмотренных типов 
пьезокерамических материалов. 

Из сравнения рис. 3, 4, аи рис. 3, 4, 6 можно заключить, что учёт неоднородности поля- 
ризации в окрестности пор обеспечивает большее убывание пьезомодулей, причём этот эффект 


сильнее проявляется для пьезомодуля ЯН: (р) . Эти зависимости лучше согласуются с известными 


экспериментальными данными [1—6] и в связи с этим можно сделать вывод о важности учёта ги- 
потезы о неоднородности поляризации пьезокерамических материалов в окрестности пор. 
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ЕРРЕСТТУЕ МОБОЧ $ТМОТАТТОМ РОВ УАВТОЦ$ ТУРЕ$ ОЕ РОВОЦ$ 
РТЕРОСЕКАМТС МАТЕВТАЕ 5$! 


А. У. Мазе4Кт 

(Роп Зае Теспткса! Упмег$Ку), 

М. $. ЗНемё5оуа 

(бощлег 5аепИЯс Сегцег, КизЗап Асадету оЁ 54епсе5) 


А сотр 5еЁ оЁЕ теепа/ соп${апЁ5 Гог рогоиз риегосегатс таепта/5 ин а/оигапсе оп пе тротодепеоих ро- 
/апгайоп п пе рогез истйу 5 а етттеа Бу те ятшайоп есптдиез юг 1е гергезетайуе уоитез, ейесвуе 
тоди/, апа Епйе еетепез. Рог йе питепса! ЧЕегтттпаНоп оЁ {те ейесвуе тоди/, 5еЁ5 о {те р/егоеесё/с чайс 
ргоМетбс Гог гергезетщайуе уоитез илён Боипаагу сопЙ#оп$ ргоуата соп&апЕ уа/иез оЁ Ве $655, зат, еес- 
{ис Не/д5, апа {те еесёис @Фэр!асетепЕ о! а Потодепеоиз Воду аге о№ме4 т {пе АМ5У$ ГЕ раскаде. Тве сасша- 
Воп5 аге ретогтеа изтд те 5реса/у аеуеюреа сотриг ргодгатз т АРОЁЕ АМ$У$ тасго [апдиаде. Тре ейесЕ 
ОР Те сей5 питьег апа иатоиз тоде[5 оЕ {те гергезещавйуе уоитез оп пе сгис/а! иа!иез о! {Те ейесйуе тойшез 
Юг а питьег о! Ёуре$ оЁ рогоиз р/егое!есё/с сегати/сс [5 пуезИдае4. ТПе гези!з соггеайоп о! {те сотршавопа! 
апа №и/-сае ехрейтепес [5 @сиеа, 1 15 зпоит {ПаЕ пе аерепаЕепсе о! те ейеснуе тоди! оп рогозйу иеаКу 
ерепсб оп {те ёуре оЁ рогоиз р/егое/есв/с таепай. 

Кеуигога5: р/егое/есёсйу, еесвоея$йсйу, рогоиз р/егоеесв"с сегат/с$, протодепеоиз роапгайоп, ейесйуе 
тоади!/, тодеб о! гергезещануе уо/итез, Ипйе-еетепЕ тетод. 





1 Тне гееагсн {5 допе мин {Пе раба! Япапса! зиррок гот Апа!уйса! БгапсН Кагдее ргодгапт оЁ {Пе Митегу оЁ ЕдисаНоп апа 
баепсе ог {Пе Виз$!ап РедегаНоп “Реувортепе ог зЧепИЙЯс роепна! оЁ МюПег $сНоо!” (ргодесЕ 1.6052.2011) апа гот ВРЕТ (12- 
07-00639_а). 
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